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ABSTRACT 

We have examined simultaneous X-ray and optical light curves of a sample of eight nearby 
Seyfert 1 galaxies observed using the EPIC X-ray cameras and Optical Monitor on board 
XMM-Newton. The observations span ~ 1 day and revealed optical variability in four of the 
eight objects studied. In all cases, the X-ray variability amplitude exceeded that of the opti- 
cal both in fractional and absolute luminosity terms. No clearly significant correlations were 
detected between wavebands using cross correlation analysis. We conclude that, in three of 
the four objects in which optical variability was detected, reprocessing mechanisms between 
wavebands do not dominate either the optical or X-ray variability on the time-scales probed. 
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1 INTRODUCTION 

Active Galactic Nuclei (AGN) emit radiation over a broad range of 
wavelengths from radio through X-rays and are variable over this 
entire spectral range. It is thought that the variability in different 
bands may be connected or driven by a single mechanism. In this 
paper we examine the potential connection between optical and X- 
ray variability in Seyfert 1 galaxies. 

Optical and UV continuum emissiorQ in AGN is believed to 
be thermal emission originating in the accretion disc surrounding 
the central supermassive black hole (SMBH). The X-ray emission, 
which shows the most rapid variability, is thought to come from 
the hot 'corona' close to the SMBH. There are two main theories 
regarding a connection between these two bands: reprocessing of 
X-rays into thermal optical emission and Compton upscattering of 
optical photons to X-ray energies. In the X-ray reprocessing sce- 
nario, X-rays from the corona are absorbed in the disc, heating it, 
which then re-radiates this energy as optical photons. The converse 
theory is that photons from the disc enter the corona and are Comp- 
ton upscattered by hot electrons to be emitted as X-ray photons. It is 
currently unclear how important these mechanisms are in forming 
the continuum emission from AGN. 

If either of these processes are important in AGN then their 
effect should be observable in the variability of the source. If op- 
tical photons drive the X-ray photons then an increase in flux in 
the optical light curve should be repeated in the X-ray band a short 
while later and vice versa. By cross correlating X-ray and optical 
light curves to look for such a correlation, we can infer if and how 
the different emission regions are connected. The traditional way 
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1 Throughout this paper we use the term 'optical' to refer to the optical-UV 
bandpass of the Optical Monitor 



to achieve this is to schedule simultaneous observations of a source 
on multiple telescopes. XMM-Newton provides a simpler way of 
performing these observations as it carries the Optical Monitor as 
well as three X-ray telescopes and can therefore perform X-ray and 
optical observations simultaneously. Here we exploit this capabil- 
ity to examine a small sample of objects and study the relationship 
between X-ray and optical variability. 

Previous studies of this subject have provided a confusing ar- 
ray of results. In some cases, no correlation is detected on short 
time-scales but light curves are correlated when averaged over a 
longer period (e.g. NGC 4051; Peterson et al. 2000). NGC 4051 
was also studied by Shemmer et al. (2003) who found evidence that 
the X-rays were led by the optical emission by ~ 2 days. For Ark 
564, on the other hand, Shemmer et al. (2001) found that the UV 
lags the X-ray emission by 0.4 days and optical lags UV by ~ 2 
days. There are also cases where the X-ray and optical light curves 
are correlated with no lag. In a 6 year observation of NGC 5548 by 
Uttley et al. (2003), a strong correlation between X-ray and optical 
light curves was detected at lag ± 15 days from data averaged at 
30 day intervals. For further examples see Uttley (2005). It is ob- 
vious that further work needs to be done to resolve the confusion 
surrounding reprocessing. 

Only two previous studies have used the simultaneous observ- 
ing capabilities of XMM-Newton to examine variability correlations 
in AGN. Mason et al. (2002) used a 1.5 day long XMM-Newton 
observation of NGC 405 1 supplemented with data from RXTE and 
found three prominent features in which the X-ray leads the optical 
emission. Using Monte Carlo simulations they derived a confidence 
of 85 per cent in these features. Arevalo et al. (2005) analysed the 
long (300 ks) observation of MCG— 6-30-15 and reported a signif- 
icant correlation with the UV leading the X-rays by 1.8 days. In 
this paper we re-examine these data along with data from six other 
Seyferts. 
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Table 1. Observation summary 



Object 


Redshift 


On time 


Revolution 


OM filter 


No. OM exposures 


P 


OM rms 


X-ray rms 






(s) 










% 


% 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1H 0707-495 


0.041 1 


46018 


0159 


UVW2 


15 


0.067 


2.0 


40.1 


Ark 120 


0.0322 


112130 


0679 


UVW2 


80 


0.051 


0.7 


2.0 


Ark 564 


0.0246 


101774 


0930 


UVW1 


49 


0.088 


0.3 


22.8 


MCG-6-30-15 


0.0077 


349286 


0301-0303 


U 


275 


<0.01 


1.9 


28.2 


Mrk 766 


0.0129 


129906 


0265 


UVW1 


75 


0.95 





25.4 


NGC 3783 (1) 


0.0097 


40412 


0193 


UVW2 


30 


0.36 


0.4 


5.8 


NGC 3783 (2) 




275633 


0371-0372 


UVW2 


53 


<0.01 


2.9 


20.7 


NGC 4051 (1) 


0.0023 


121958 


0263 


UVW1 


43 


<0.01 


1.5 


47.6 


NGC 4051 (2) 




51866 


0541 


UVW1 


15 


<0.01 


0.8 


43.8 


NGC 7469 


0.0163 


164128 


0912-0913 


UVW2 


111 


<0.01 


1.9 


7.7 



The remainder of this paper is organised as follows. In sec- 
tion|2]we describe the observations selected and the data reduction 
procedures. The robustness of our technique is discussed in sec- 
tion[3] Section[4]contains the results of the analysis and in section[5] 
we compare our results to those made previously and discuss the 
implications of these results to the study of AGN variability. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Sample definition and observations 

The satellite XMM-Newton carries three X-ray telescopes and a co- 
aligned Optical Monitor (OM; Mason et al. 2001), comprising a 
30 cm Ritchey-Chretien telescope with a micro-channel plate in- 
tensified CCD. This gives a field of view of 17 arcmin square and 
spatial resolution of ~ 1 arcsec, with 0.48 arcsec pixels. The OM 
can take images using any of six narrow-band filters sensitive over 
the range 1800 — 6000 A, thereby providing optical imaging and 
timing information simultaneous with the X-ray observations. The 
combination of long X-ray exposures with simultaneous optical 
data makes XMM-Newton ideal for probing the relation between 
X-ray and optical variability on short time-scales. 

The XMM-Newton Science Archive (XSA) contains all avail- 
able data products from the satellite. The XSA was used to find 
long observations (> 30 ks) of nearby (z < 0.1) Seyfert 1 galax- 
ies. Objects were selected if there were 15 or more exposures taken 
with a single filter of the OM. Eight objects were found to satisfy 
these criteria: 1H 0707-495, Ark 120, Ark 564, MCG-6-30-15, 
Mrk 766, NGC 3783, NGC 4051 and NGC 7469. Four of these 
sources had observations in two or more revolutions in the same 
filter. A summary of the observations of these sources is given in 
Table Q] 

2.2 Predicted time delays for the sample 

The eight targets selected all have published black hole mass esti- 
mates, which allowed us to make simple predictions for the inter- 
band time delays that might be expected from reprocessing models. 
Firstly we assumed all heating in the disc is generated by the vis- 
cosity of the disc itself. Different radii of the disc are at different 
temperatures and can be modelled as blackbodies with a particular 
peak wavelength. Using the effective wavelength of the appropriate 
OM filter, we calculated the temperature of a blackbody spectrum 



that peaks in the observed band using Wien's Law. This tempera- 
ture was then used with the standard formulae for a geometrically 
thin, optically thick accretion disc to determine the approximate 
radius from which the emission originates: 

"■-" 4 (£r 

where MIM e is the accretion rate as a fraction of the Eddington 
accretion rate, Mg is the mass of the black hole in units of 10 8 Mq 
and r/Rs is the radius in terms of the Schwarzchild radius (Peter- 
son 1997). Using values of 0.01 and 0.1 for M/Me, and published 
black hole masses (see Table [2}, we derived a range for the radius 
of the optical region being observed. Assuming the X-ray emission 
comes from close to the central black hole, the light travel time to 
the radius of optical emission region gives a rough estimate for the 
time delay between these two regions. 

As another check on the delay expected, we used the assump- 
tion that the disc is heated entirely externally by X-rays. In this 
case, the energy of the blackbody peaking in the OM filter is com- 
parable to the total X-ray luminosity. Using the temperature as cal- 
culated above, Stefan's Law was used to determine the radius of the 
optical emission region: 

L x = 4-Kr 2 aT 4 (2) 

where a is Stefan's constant. The X-ray luminosity, Lx, is obtained 
from the 0.1 — 10 keV X-ray spectrum. Again the radius was con- 
verted to an approximate time delay in light days. The results of 
these calculations are shown in Table [2] The majority of sources 
are expected to show time delays of less than a day which should 
theoretically be observable with XMM-Newton. 



2.3 Optical light curves 

The Observation Data Files (ODFs) for each source were extracted 
from the XSA and processed using the task omichain of the 
XMM-Newton Science Analysis System (SAS v6.5.0). Omichain 
takes the ODFs and performs all standard OM processing steps on 
them to produce the final products - images and source lists. Once 
the images have been created, omichain performs source detec- 
tion and photometry routines. Finally the source positions are con- 
verted from detector to sky coordinates and a combined source list 
from all exposures is created. 
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Table 2. Estimated time delays 



Object 




Int. delay 




Ext. delay ] 




(10 6 M ) 


(days) 


(10 44 ergs" 1 ) 


(days) ' 


(1) 


(2) 


(3) 


(4) 


(5) ; 


1H 0707-495 


2.3 1 


0.02-0.03 


0.42 


i 

0.16 


Ark 120 


150±19 2 


0.25-0.54 


20.1 


1.10 ' 


Ark 564 


1.09 1 


0.01-0.03 


4.5 


0.98 : 


MCG-6-30-15 


9 n+1.83 


0.03-0.07 


1.14 


0.69 ; 


Mrk 766 


0.63 4 


0.01-0.02 


1.83 


0.62 


NGC 3783 


29.8±5.4 2 


0.09-0.18 


2.55 


0.39 l 


NGC4051 


1.91±0.78 2 


0.02-0.04 


0.03 


0.08 « 


NGC 7469 


12.2±1.4 2 


0.05-0.10 


3.9 


0.48 1 



Peterson et al. (2004); 3 McHardy et al. (2005); 



1 Zhou & Wang (2005); 
4 Botte et al (2005) 

NOTES: Values in column 3 are based on the assumption of internal disc 
heating by viscosity. The lower limits use a value of MIMe = 0.01 and the 
upper limits use 0.1. X-ray luminosity (column 4) is for the band 0.1-10 keV. 
Values in column 5 are based on the assumption of external disc heating by 
X-rays. 



To create an optical light curve, photometry must be per- 
formed on each exposure of each object. Although omichain au- 
tomatically performs photometry on each image, the presence of 
host galaxy emission in some images causes the routine to flip be- 
tween 'extended' and 'point source' modes of photometry, and so 
does not perform a uniform extraction of all the images. As the 
aim is to obtain accurate photometry of the nucleus and not include 
the host emission this is obviously undesirable. In order to ensure 
uniform data reduction we used custom made scripts to perform 
aperture photometry and apply the appropriate OM corrections to 
the resulting data, modelled closely on the SAS processing tasks 
omdetect and ommag. 

Simple aperture photometry was applied to the images pro- 
duced by omichain (after the modulo-8 correction has been ap- 
plied by the task ommodmap). Based on the user defined source 
position (in detector coordinates), the exact source location was 
obtained by repeatedly centroiding on the counts within a small 
box, then the counts within a circular aperture were accumulated. 
For all sources, except MCG— 6-30-15, a circle of radius 12 pix- 
els was used as the source aperture to match the calibration (see 
below). In the case of MCG— 6-30-15, a 12 pixel radius aperture 
overlapped a nearby star so, to avoid this, an aperture of 7 pixels 
was used for all exposures of MCG— 6-30-15. For most sources, an 
annulus placed around the source aperture was used to extract the 
background count rate. The background annulus was placed at least 
30 pixels out from the source to avoid the extended PSF of the UV 
filters. In the cases of Ark 564, MCG-6-30-15 and NGC 4051 the 
host galaxy is visible around the AGN. As the host galaxy should 
not be varying in flux and it is the flux variability, not absolute flux, 
that we are studying, we chose to attempt no subtraction of the host 
galaxy flux from the flux of the central AGN. To get an accurate 
value for the actual background of the field, a circular background 
aperture away from the host galaxy was used in these cases. 

In order to produce accurate count rates from aperture pho- 
tometry of OM images there are five corrections that must be 
applied - point spread function (PSF1), coincidence loss (CL), 
deadtime (DT), a UV specific PSF correction (PSF2) and time- 
dependent sensitivity degradation (TDS). Coincidence loss is a loss 
of flux that occurs when more than one event is collected in the 



same pixel in a single CCD frame, and the CL correction is a non- 
linear correction that has been calibrated empirically and found to 
be self-consistent when applied to an aperture radius of 6 arcsec (12 
pixels). Before the CL correction can be applied, the PSF1 correc- 
tion must therefore be implemented to scale the count rate from the 
aperture used to an effective aperture radius of 12 pixels. When the 
CCD is read out there is a short amount of time, the frame transfer 
time, when any incoming photons are not detected. The deadtime 
factor corrects for this loss. In the UV filters the wings of the PSF 
are more extended than in the optical filters so a further correction 
is made to count rates from exposures in these filters to account for 
the loss of count rate from the wings. Finally the time-dependent 
degradation factor takes account of the fact that the performance of 
the OM gradually decreases over time. The uncertainties on each 
flux point were derived by propagating the Poissonian error on the 
source and background count rates. 



2.4 X-ray light curves 

X-ray light curves were produced using data from the European 
Photon Imaging Camera pn camera (EPIC-pn) onboard XMM- 
Newton. Event lists for each observation were extracted using the 
standard processing procedure, epchain. In order to check for 
background flaring, 10 — 15 keV light curves were extracted from 
off-source regions and examined by eye for periods of high, vari- 
able background. Where such intervals occurred they were re- 
moved by filtering out all time intervals during which the count 
rate in the high energy light curve exceeded a critical value that 
was chosen to be slightly above the quiescent level, but was 
slightly different for each observation. After filtering out back- 
ground flares, source light curves were extracted from a 40 arc- 
sec circular region using only single or double-pixel events in the 
range 0.2 — 10.0 keV, and background subtracted using a nearby 
source-free region. The light curves were binned such that the X- 
ray sampling is comparable to the sampling of the OM light curves. 

Figure Q] shows the optical/UV and X-ray light curve for the 
sources with detectable X-ray and optical variability. MCG— 6-30- 
15, NGC 3783 and NGC 7469 all had consecutive observations 
which have been combined. 



3 TESTING THE METHOD 
3.1 Monte Carlo tests 

In order to test the robustness of our light curve extraction (out- 
lined above) to PSF losses and host galaxy contamination, random 
data were simulated and subjected to the same aperture photometry 
procedure. As the spacecraft pointing can change during an obser- 
vation by as much as a few arcsec, causing the target source to drift 
several pixels across the image, this effect was also included in the 
simulations. The simulations were designed to match the MCG— 6- 
30-15 data, chosen because it is the longest dataset and shows the 
most contamination from its host galaxy. The simulation procedure 
was as follows. 

The 'true' source position on a 224 x 224 pixel detector grid 
was defined as (xq, yo) based on the average detector position of 
MCG— 6-30-15 in the real OM images. The effect of random point- 
ing changes was simulated by giving a random offset to the precise 
position of the point source in each image. A point source was then 
added to the image at the new position (xo + Ax, yo + Ay) using 
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Figure 1. X-ray (0.2 - 10 keV) and optical light curves of MCG-6-30-15 (Rev. 0301-0303), NGC 3783 (Rev. 0371-0372), NGC 4051 (Rev. 0263), NGC 
4051 (Rev. 0541) and NGC 7469 (Rev. 0912-0913). X-ray light curves are binned to match the optical sampling. Error bars are included for all count rate 
measurements, but for several of the X-ray light curves the error bars are smaller than the data point symbols. 



the PSF model obtained from the latest calibration files. In prac- 
tice a differential PSF model was calculated from the cumulative 
PSF model listed in the calibration files, and linearly interpolated 
to produce a function sampled with 0.1 pixel resolution out to 12 
pixels. This was used to define the image of the point source from a 
smooth PSF projected onto a relatively coarse pixel grid. The host 
galaxy emission was modelled as an ellipsoid with an exponen- 
tial surface density profile, with spatial parameters and normalisa- 
tion (relative to the point source) taken from the spatial analysis of 
Arevalo et al. (2005), kindly supplied by Dr Arevalo. The source 
image was then normalised to match the typical counts per image 
from the observations, and a typical (flat) background level (~ 3.3 
ct/pix) was added. The counts per pixel were then randomly drawn 
from a Poisson distribution to simulate random counting statistics. 
The simulated images were then subjected to aperture photometry. 
Specifically, repeated centroiding was used to find the exact source 
position in each image, then counts were accumulated from a cir- 
cular region and the PSF1 correction was applied exactly as for 
the real data. The effects of the other corrections, namely CL, TDS 
and DT, were not included. The latter two corrections are simply 
constant factors determined for each exposure, and the CL is well 
calibrated for aperture photometry (assuming the PSF1 correlation 
is correct). 

The simulations were designed to reproduce the effects of 
aperture photometry on a constant point source, the position of 
which is slightly different for each image, superposed on a real- 
istic host galaxy image and flat constant background level, all sub- 
jected to Poisson variations. The end result of 200 simulations was 
a time series that could be fitted with a constant to test for spu- 
rious variability. Using random position offsets (As, Ay) drawn 
from a Gaussian of mean zero and standard deviation 0.2 pixels, 
the result was a fit statistic of x = 211.29 with 199 degrees of 
freedom (dof), a perfectly acceptable fit (p = 0.262), indicating no 
statistically significant variability was present given 200 images. A 
longer simulation using 5000 images did show significant variabil- 
ity (x 2 = 5432.56 for 4999 dof, p = 1.2 x 10" 5 ) but only at the 
0.4 per cent level, which is too small to be detected in less than 200 
images. This excess variability is also far smaller than the measured 
amplitudes in any of the observations (see below). 

Allowing for a 10 pixel drift in both x and y directions over 
the course of 200 images, comparable to the true drift during the 
MCG— 6-30-15 observation, plus a random offset for each image, 
gave similar results: x 2 = 205.19 for 199 dof (p = 0.39). There 
was no statistically significant correlation (r = 0.044, p — 0.53) 
between the measured flux and the offset between the centroid and 
the stationary 'first guess' position (xo,yo) given to the photometry 
algorithm, indicating that the photometry method is not biased by 
even fairly large changes in spacecraft pointing. The same results 
were obtained when the host galaxy was not included in the simu- 
lations, proving that contamination by host galaxy light within the 
source aperture does not affect the variability. Changing the ran- 
dom offsets (Ax, Ay) to have a standard deviation of 0.5 pixels in 
each direction also gave entirely consistent results. 

3.2 Empirical tests using field stars 

Satisfied that our method performed well using simulated data we 
then tested it on actual data from the OM. For several of the sources 
in the sample (Ark 564, Mrk 766, MCG-6-30-15), stars were vis- 
ible in the field of view. To prove that the method does not intro- 
duce spurious variability, the process was tested on a selection of 
these field stars. The field around Ark 564 contains many objects. 
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Table 3. Field sources used for testing 



Field target 


RA 


Dec 


V 




(degrees) 


(degrees) 




(1) 


(2) 


(3) 


(4) 


Ark 564 


340.593 


29.725 


0.02 




340.701 


29.756 


0.70 


MCG-6-30-15 


203.966 


-34.137 


0.21 


Mrk 766 


184.611 


29.813 


0.80 


NGC4051 


180.786 


44.539 


0.21 




180.793 


44.541 


0.56 




180.763 


44.541 


0.77 



For our field star tests we chose two objects of similar brightness 
to Ark 564 which were far enough removed from other objects in 
the field to allow accurate light curves to be extracted. The other 
fields examined here contained few objects other than the target so 
all field stars (unless contaminated by being too close to the target) 
were tested. 

Table [5] summarises the results of these tests. A p value of 
^> 0.01 indicates that a light curve is non-variable. Three of the 
four stars tested proved to be non- variable proving that the method 
used here does not introduce spurious variability. One of the stars 
tested in the Ark 564 field has a p value of 0.02 implying possible 
variability in the source. We performed CCF analysis between this 
object and the Ark 564 light curve to ensure the detected variability 
is not a systematic problem also affecting Ark 564. No correlations 
were found between the two objects. We suggest that the star is in- 
trinsically variable and does not cause a problem with our method. 

As another test of whether or not the method is robust, light 
curves for extended sources in the NGC 4051 host galaxy were 
extracted. These will not vary on the time-scales probed here so 
a non- varying light curve should definitely be produced. The light 
curves extracted did indeed prove to be non- varying confirming that 
the method does not introduce spurious variability. 



4 RESULTS 

To determine whether or not the optical light curve is variable, a 
X 2 test was performed on each light curve. The resulting p values 
are given in Table [T] Values of p < 0.01 indicate variation at the 
99 per cent confidence level or greater and values between 0.01 
and 0.05 are possible detections of variability (95 — 99 per cent 
confidence). Variability was detected in the OM light curves at 99 
per cent confidence in four out of eight objects. 

Also shown in Table Q] is the optical and X-ray rms variabil- 
ity amplitude for each source. The typical error in these values is 
around one per cent so, depending on the brightness of the source, 
an object with rms amplitude greater than one should be detectable 
as varying. The optical rms values shown here are probably un- 
derestimated as they include the total flux of the source with no 
host galaxy subtraction. For the four sources that show significant 
variability in both optical and X-ray bands, the rms variability am- 
plitude was calculated in luminosity units, as tabulated in Table [4] 

A standard tool for testing and measuring correlations be- 
tween two time series is the cross correlation function (CCF), 
which gives the degree of correlation between the two series as 
one series is shifted backwards and forwards in time relative to the 
other. Unfortunately, when the data are not exactly evenly sampled 
one cannot apply the standard, 'textbook' formulae for estimating 
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Table 4. RMS variability 



Object 


OM rms 


X-ray rms 




(10 42 ergs" 1 ) 


(lO^ergs- 1 ) 


(1) 


(2) 


(3) 


MCG-6-30-15 


0.2 


35 


NGC 3783 (2) 


3.3 


51 


NGC4051 (1) 


0.08 


2.9 


NGC 4051 (2) 


0.5 


0.7 


NGC 7469 


7.3 


30 



the CCF (Box & Jenkins 1976; Priestley 1981; Brockwell & Davis 

1996) . There are two reasons why the present data are not evenly 
sampled: the OM images are rather evenly spaced in time but there 
are a few gaps, and observations spanning more than one XMM- 
Newton revolution will be split by the orbital gap. Suggested solu- 
tions to the problem of uneven sampling include interpolating the 
data onto an even grid, or calculating the correlation using pairs 
of data points that have a lag falling within a finite range (i.e. 
Tij = t x (i) — t y (j) falls in the interval [r, r + At]). Specific algo- 
rithms that have been used in AGN monitoring include the interpo- 
lated cross correlation function (ICCF; Gaskell & Peterson 1987), 
the discrete correlation function (DCF; Edelson & Krolik 1988) and 
the z-transformed discrete correlation function (ZDCF; Alexander 

1997) . All three of these were used in the present analysis and the 
results were found to be consistent. 

A second problem that plagues CCF estimation from AGN 
data is that the data are often non-stationary, meaning that the mean 
and/or variance of a short light curve will change over time, due to 
variations on time-scales longer than the length of the observation. 
Following the advice of Welsh (1999) two precautions were used to 
mitigate the side-effects of non-stationarity. Firstly a linear function 
was fitted to the data (using ordinary least-squares) and removed to 
suppress any long time-scale variation, and secondly the mean and 
variance were calculated 'locally', using only the data contributing 
to a particular lag. Following another suggestion in Welsh (1999) 
the correlation coefficient was computed using both the raw data 
and the ranked data (each data point being correlated was replaced 
with its rank). In principal using the ranked data should improve 
sensitivity to non-linear correlations and make the statistic more 
robust to outliers, compared to using the raw data. In practice the 
two were found to be in close agreement. 

Confidence limits on the CCF were computed using Bartlett's 
formula (Bartlett 1955) as an approximate but simple test of the 
significance of any correlations. In particular, standard deviation 
of the distribution of CCF points around an expected value of zero 
was computed using equation 1 1 . 1 .9 of Box & Jenkins ( 1 976). This 
assumes that in reality the two series are not correlated, and can be 
used to test against this hypothesis. Under the assumption of no 
cross-correlation (i.e. the null hypothesis that the expectation value 
of the CCF is zero) Bartlett's formula may be written as: 

cr 2 {CCF(rj)} = -jj— VACF 1 (r i )ACF 2 (r < ) (3) 

i 

where the two ACF functions are the empirical auto- 
correlation functions and the summation is carried out over all mea- 
sured lags. Here Apair is the number of pairs of data points contri- 
bution at the lag Tj . 

The 95 per cent confidence region was calculated as ±1.96cr. 



The confidence region for the correlation coefficient at a given lag 

— 1/2 

scales with N ^ . The bounds are smallest at small time delays, 
where the two time series overlap the most, and increase towards 
the edges of the CCF. The CCFs were examined over a range of 
lags spanning, in either direction, half the observation duration. At 
longer lags the number of pairs that contribute to the CCF calcula- 
tion rapidly decreases, leading to highly uncertain CCF estimates. 
It should be noted however that Bartlett's formula is strictly only 
valid when the number of data points N is large and the data are 
stationary. As the X-ray/optical data do not satisfy either condi- 
tion well, the confidence limits should be considered as rough es- 
timates of the size of the true confidence region. Bartlett's formula 
is discussed in most standard references for time series analysis in- 
cluding: Jenkins & Watts (1968; sect 8.2.1), Box & Jenkins (1976; 
sect 11.1.3), Priestley (1981; sect 9.5), Bendat & Piersol (1986; 
sect 8.4), Kendall & Ord (1990; sect 11.3), or Brockwell & Davis 
(1996; sect 7.3.4). 

Previous papers discussing inter-band correlations in AGN 
variability have made use of Monte Carlo simulations to test the 
significance of any possible correlations. In principle the Monte 
Carlo approach should be robust but in practice it does rely on the 
power spectra of the two light curves being well understood in or- 
der to produce realistic simulations. Technically the CCF is not a 
'pivotal' statistic, meaning it is sensitive to 'uninteresting' param- 
eters such as the power spectral shape, and so simulations must be 
made as carefully as possible. For the present analysis it was not 
possible to obtain reliable power spectra from the limited OM data, 
and so Bartlett's formula (which uses the empirical auto-correlation 
function) was used as a guide of significance. It is also worth noting 
the results of Monte Carlo tests are sensitive to which of the pair of 
light curves are randomised - whether the X-ray, or optical, or both 
dataset are randomised can substantially affect the apparent signifi- 
cance, and as yet there is no accepted preference. This choice is not 
required for Bartlett's formula to be used. 

Figure [2] shows the CCF plots for the objects studied which 
displayed optical and X-ray variability. A positive lag in the plots 
implies the X-ray variations lead those in the optical and vice versa 
for a negative lag. The 95 and 99 per cent confidence limits are in- 
dicated by dotted and dot-dashed lines respectively. In the case of 
NGC 4051 (Rev. 263) and NGC 3783 the 95 per cent confidence 
region is sufficiently broad that these data do not impose useful 
constraints on the possibility of a correlation. The error bounds, as 
computed using the Bartlett formula, were large in these cases be- 
cause the empirical ACFs for both optical and X-ray time series 
were very broad. Particularly in the case of NGC 3783 the X-ray 
power spectrum over the time-scales studied here is rather steep, 
leading to smooth light curves and a broad ACF, and hence a large 
error in the Bartlett formula. This accounts for the fact that real- 
isations of processes with steeper power spectra (or equivalently, 
broader ACFs) are more 'weakly non-stationary'. However, in the 
cases of NGC 4051 (Rev. 541), NGC 7469 and MCG-6-30-15 the 
CCF is sufficently well constrained to exclude a strong correlation, 
i.e. 95 per cent confidence limit of | CCF(r) | < 0.5, over a plausi- 
ble range of r. 

Some previous studies (e.g. Nandra et al. 2000) have detected 
a correlation between the optical light curve and X-ray photon in- 
dex but no correlation between optical and X-ray light curves. In 
light of this we also cross correlated the optical light curves of the 
objects studied here with the 0.7 — 2/2 — 10 keV hardness ratio 
(which should track the photon index) but found the cross corre- 
lation functions to be very similar to the ones between X-ray and 
optical flux. 
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Figure 2. CCF plots of MCG-6-30-15 (Rev. 0301-0303) NGC 3783 (Rev. 0371-0372), NGC 4051 (Rev. 0263 - 95 per cent confidence limit just visible), 
NGC 4051 (Rev. 0541) and NGC 7469 (Rev. 0912-0913). The DCF is plotted in black and the ICCF is the red dashed line. The dotted and dot-dashed lines 
represent the 95 and 99 per cent confidence limits respectively. 



5 DISCUSSION 



5.1 Summary of results 

Using data from the OM and X-ray telescopes aboard XMM- 
Newton we have examined a sample of eight objects to look for ev- 
idence of reprocessing between wavebands. Optical variability on 
short time-scales was detected at 99 per cent confidence in four out 
of eight objects and rms variability amplitudes in all cases show that 
the X-ray variability is stronger, in both fractional terms and abso- 
lute luminosity units, than the optical variability. In three out of the 
four objects showing detectable optical variability, the cross corre- 
lation function is sufficiently well constrained to exclude a strong 
X-ray/optical correlation and hence reprocessing on time-scales of 
< 1 day is not a dominant mechanism. 



5.2 Comparison with previous studies 

5.2.1 Observations using XMM-Newton data 

The observation of NGC 405 1 in Revolution 0263 was also inves- 
tigated by Mason et al. (2002) for inter-band correlations. They 
supplemented the XMM-Newton data with X-ray observations from 
RXTE and their CCF analysis covers positive lags up to two days. 
Over the region that their CCF plot overlaps with ours (0—0.6 days) 
the results are consistent. We observe a peak in the CCF around 0.2 
days with a value of 0.75 similar to their peak of ~ 0.6. Using 
simulations of uncorrelated red-noise light curves they place a con- 
fidence of 85 per cent in the correlation. Our peak lies between 
68.3 and 95 per cent confidence levels, certainly not significant at 
99 per cent. However, as mentioned in section|4] the confidence re- 
gions in this case are not well enough constrained to comment on 
the correlation. 
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The OM observation of MCG— 6-30-15 was previously anal- 
ysed by Arevalo et al. (2005), who reported a positive correlation 
with DCF max = 0.82 at 160 ks (around 1.8 days) with the UV 
leading the X-ray data. Using Monte Carlo simulations they esti- 
mated a confidence limit of 98.5 per cent (i.e. a p-value of 0.015). 
By contrast our analysis found a much lower CCF peak (0.35) 
which was not considered a significant detection of correlation. 

There are several differences between the analysis of Arevalo 
et al. (2005) and that of the present paper. In terms of the light curve 
extraction, Arevalo et al. (2005) performed photometry on attitude 
corrected images by fitting an empirical PSF model to the source 
image, while simultaneously modelling the host galaxy emission. 
The present analysis used aperture photometry on (detector space) 
images based on the SAS processing pipeline tasks, including cor- 
rections for coincidence loss, deadtime and the counts/frame de- 
pendence of the PSF model, none of which were included by 
Arevalo et al. (2005). These differences in extraction technique 
explain the apparent difference between the two reductions of the 
same data. Arevalo et al. (2005) also used coarser binning of the 
light curves, but applying the same binning to our light curves did 
not substantially decrease the discrepancy. As shown in Section [3] 
the aperture photometry method employed in this paper is robust to 
host galaxy presence and to pointing changes but, unlike the anal- 
ysis of Arevalo et al. (2005), our method includes all the known 
time-dependent flux corrections that should be applied to OM data. 
Therefore, in principle at least, this reduction should be the more 
reliable. 

In practice most of the difference between the correlation 
strengths is due to the more conservative approach employed in the 
present analysis for calculating the CCF. Arevalo et al. (2005) pre- 
sented the DCF after normalising by the noise-subtracted variances 
(i.e. the variances used in the denominator of the correlation coef- 
ficient formula had the mean square error subtracted). The present 
analysis used ICCF, DCF and ZDCF methods without the noise- 
subtraction term, employing linear and rank-order correlation coef- 
ficients, the latter should be more robust to outliers. Including the 
noise-subtraction term increased the CCF value by « 0.2 at the 
claimed lag, but this only affects the absolute value of the CCF and 
not its significance (since the effect should be calibrated out during 
the significance test). 

A second difference is that, in the present analysis the time se- 
ries were 'detrended' before calculating the CCF by subtracting off 
a linear function. As recommended by Welsh (1999) this helps re- 
move any bias caused by very long time-scale variations manifested 
as quasi-linear trends over the data that can lead to spuriously high 
correlation coefficients. Applying the detrending reduced the CCF 
by ~ 0.2 at the claimed lag. If the correlation at ~ 1.8-d was robust 
it should not be affected by removal of longer time-scale (~ 5-d) 
variations. 

J. 2. 2 Observations from other telescopes 

Several other sources discussed here have also been studied previ- 
ously: Ark 564, NGC 4051 and NGC 7469. 

Ark 564 was observed by Papadakis et al. (2000) over a period 
of 21 days at the Skinakas Observatory in optical and simultane- 
ously in X-rays with ASCA. Papadakis et al. (2000) probed slightly 
longer time-scales than we do but, in line with the trend reported 
here, they find no significant correlations between the two bands. 
Shemmer et al. (2001) examined Ark 564 over a longer period (50 
days) in UV with HST and X-ray with RXTE and ASCA . They found 
evidence of UV emission following X-ray by ~ 0.4 days. 



Shemmer et al. (2003) studied NGC 4051 and found evidence 
for a correlation between X-ray and optical bands with a lag of less 
than one day. Their campaign covered a period of three months 
with approximately daily observations and so probed a longer time- 
scale than we examine here. NGC 4051 was also investigated by 
Peterson et al. (2000). On short time-scales they found no evidence 
for correlations. Having smoothed the data in 30 day bins however 
a correlation consistent with zero lag is detected. 

A 30 day observation of NGC 7469 by Nandra et al. (1998) 
found that at least two variability mechanisms are required to de- 
scribe their result. Cross correlation of UV and X-ray light curves 
gave evidence for a lag of around four days. This lag explained 
the differing times of maxima in the two light curves but failed 
to describe how the minima in the light curves occurred simulata- 
neously. 

As none of the studies described above probe the same short 
time-scales examined here it is difficult to make a comparison be- 
tween them and the results described in this paper. An overall pic- 
ture seems to be emerging of no correlation between X-ray and 
optical emission over sub-day periods but significant correlations 
only appearing on much longer time-scales. 

5.3 Implications for AGN variability 

In line with other work we find that X-ray rms variability ex- 
ceeds that in the optical in all cases. However, in the objects with 
well constrained CCF (MCG-6-30-15, NGC 4051 (Rev. 0541) and 
NGC 7469), the absence of a strong correlation implies the major- 
ity of the variance in the optical band on short time-scales is un- 
correlated with the short time-scale X-ray variations. Papadakis et 
al. (2000) suggest that this behaviour could be explained if only 
a small section of the accretion disc is visible to the X-rays or if 
the X-ray emission is anisotropic. Two of the observations show- 
ing optical variability, NGC 3783 and NGC 4051 (Rev. 0263) have 
poorly constrained CCF confidence regions and so a firm conclu- 
sion on these objects can not be drawn. 

In terms of the overall picture emerging from variability stud- 
ies this solidifies earlier findings (Edelson et al. 2000; Uttley et 
al. 2003) that, on short time-scales, X-ray and optical variations in 
Seyfert 1 galaxies are uncorrelated. Many papers (see section [5T2l > 
find some evidence for correlations on longer time-scales. It seems 
that rapid X-ray and optical variations in Seyfert 1 galaxies occur 
through separate mechanisms but long-term variations are some- 
how connected. 

Intensive observations with XMM-Newton are not suitable for 
probing longer time-scale variations. However, one promising al- 
ternative is the Swift Gamma Ray Burst Explorer. Swift carries co- 
aligned X-ray (XRT; Burrows et al. 2005) and UV/optical (UVOT; 
Roming et al. 2005) telescopes covering essentially the same band- 
pass as XMM-Newton, but is in a low-Earth orbit making it bet- 
ter suited to regular 'snapshot' monitoring of sources over weeks- 
months. The combination of XMM-Newton and Swift is therefore 
capable of probing X-ray/optical correlations with unprecedented 
detail on all time-scales from hours to weeks. 
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